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We have given a straightforward method to solve the problem of noncentral an- 
harmonic oscillator in three dimensions. The relative propagator is presented by 
means of path integrals in spherical coordinates. By making an adequate change of 
time we were able to separate the angular motion from the radial one. The relative 
propagator is then exactly calculated. The energy spectrum and the corresponding 
wave functions are obtained. 
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I. INTRODUCTION 

As we know the Feynman path integral formalism has played an important role in the 
comprehension of the quantum phenomena [IH3]. In nonrelativistic quantum mechanics, 
many problems have been exactly solved starting from their classical origins by the use of 
path integrals jH |5]. Therefore, it has been proved that this method is powerful in finding 
exact propagators and wave functions and studying quantum behavior of several systems [BJ . 
At a certain epoch, however, the development of the path integral method wasn't rapid com- 
pared with Schrodinger and Heisenberg approaches. The introduction of the space and time 
transformation to solve the nonrelativistic Coulomb problem and other basic potentials was 
the decisive stage for its development [74 IE]- Having crossed this crucial stage, the path inte- 
gral quantization method becomes at the present time more important than the Schrodinger 
equation. Its rise unveils itself from its applications in several fields of physics such as quan- 
tum mechanics [TJ-[B], statistical physics [8], quantum field theory [9j [10], condensed matter 
[IT] , cosmology [121 US] and black hole physics [TS] . We can find also a modern introduc- 
tion to path integrals with significant applications in [IB]. This explains the increasing 
interest in the development of path integration techniques. 

However to our knowledge there is no report on the path integral treatment for a class 
of noncentral potentials which are of interest in understanding the nuclear shell structure 
within the framework of a modified mean field model [17]. These potentials that describe 
the rotational- vibrational motion of the nuclear system, have the form 

V(r,9)= 1 -^ + q^ (1) 

where ^fiu 2 r 2 is the usual three dimensional harmonic oscillator that explains the magic 
numbers for spherical symmetric nucleus, u (6) is a function of 9 and q is a deformation 
parameter. We note that a particular case in two dimensional space is discussed in [18J. 
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In the present paper, we use of the Feynman path integral method to solve the problem 
of a nonrelativistic particle subjected to the following noncentral potential 

VM = -V + -^ + — 2 + ^^ + ^^, (2) 

where Vq, a, (5 and 7 are constants. We note that some particular cases of this potential 
are studied in [T9H2T] . Also some similar potentials have been recently discussed either in 
nonrelativistic quantum mechanics or in relativistic theory [22H26] . 

In the first stage we formulate the problem in spherical coordinates path integrals. Then 
by the use of a temporal transformation we separate the angular motion and we do integra- 
tion over 6 and p to obtain the exact propagator. Finally, we extract the energy spectrum 
and the wave functions. 



II. PATH INTEGRAL FORMULATION 

In nonrelativistic quantum mechanics, the propagator or the transition amplitude from 
initial state \ f a ) to final state \r b ) for a physical system governed by the Hamiltonian H and 
the Lagrangian C, is defined, in configuration space, by a matrix element of the evolution 
operator 



K (f b ,t b ;f a ,t a ) 



n 



exp 



h 



H (t b - t a ) 



which admits the following functional integral representation 



K (n, t b ; f a , t a ) = I Df(t) exp ( -S [f(t)\ ) , 



(3) 



(4) 



where the action 5* [r(t)] is a functional of continuous trajectory f(t) connecting space-time 
point (r a ,t a ) with (f b ,t b ) 



with the Lagrangian 



and the potential 



S [r(t)} 



L = ^(x 2 + y 2 + z) -V(x,y, 



dtL {f , r) , 



(5) 
(6) 



^ + ^ (7) 



V (x, y, z) = -V + \^ (x 2 + , 2 + , 2 ) + 2 ^ 2 / y l h +z2) + 2/i {x2 + y2) ■ w 
In its discrete form K (f b , t b ; r a , t a ) can be written as 



K (f b , t b ; f a , t a ) = K (f fc , f a ; T) = j Dx j Dy j Dz exp | A n 



(8) 



where 



.4, 



/1 {Ax n ) 2 + {Ay n f + {Az r , 



eV (x n , y n , z n ^j 



(9) 



2iirhe 



and 

/« . .N-l .N-l .N-l . N 

Dx Dy Dz=Vm \{dx n \{dy n \{dz n / J] 

7 7 e^O 7 71=1 7 71=1 7 71=1 J 71=1 

Here we have used the standard notation 

Ag n = g„ - g„-i 

gjv = q (h) 

Qo = q {to) 

T = t b -t a = Ne. 

Let us now search for a path integral representation in the spherical coordinates 

x n = r n sin 9 n cos y? n 
y n = r n sin 9 n sin y9„ 

Z n = T n COS $ n . 

It is obvious that the measure 

N-l .N-l .N-l N 



\\dx n I J | dy n j J J r/;„ ] [ 



71=1 71=1 71=1 71=1 

will be in spherical coordinates 

N-l .N-l „N-1 N-l 



2inhe 



f n dr n I n *>» / n ^ n (dy 3 * 

71=1 71=1 71=1 71=1 



and the propagator takes the following form 



, N-l „N-1 „N-1 N-l 



K (f b , f a , T) = Jim j H dr n j \[ d6 n f \[ d<p n J] r 2 n sin 6 r 

e^0° ^ n=l n=l ^ n=l n=l 

E( i li r 2 n + rl_ ± - 2r„r n _i cos n cos 6 n -i 
[h2 e 

% l±r n r n _ x sin 6 n sin i 
cos (A<p n ) - -eV (r„, 6> n ) 

ft £ ft 

Taking into account that 

AT-l AT 

n r « = — n ( r « r - 1 ) 

and using the following formulae 

exp (-i® cos (Aip n fj =^2l m e 



2inhe 



n=l 



+oo 

imA(p n 



where I m (x) is the modified Bessel function, we obtain 

K (f b , f a ; T) = 



r b r a Vsin 9 b sin6» a 



N-l 



. N-l N 



n=l n=l n=l n=l 

N +00 



n=l m=— 00 



exp 



i /i r„ + r„_! - 2r„r n _i cos 0„ cos 6> n _i i 



h2 



h 



V(r n ,e n )\e\ . 



For small £, the modified Bessel function behaves as follows 



Im ~i 



what permits to us to write 
K(f b ,f a ;T) 



IE 

2na 



1 



x exp 



.a is 
~ % ~E ~ 2a 



m — 



2nr b r a \J sin 9 b sin 9 a 

N-l .N-l N 



l\dr n / I I <W„ J] ( 

^n^'n— 1 ) 

n=l n=l n=l 

N +00 

n e ex p 



/i 



2inhe 



2N 



N-l 

n 

n=l 



2tt 



n=l m=—oo 



i V r l + r l-i - 2r n r„_i cos (A0 n ) 



(m 2 - |) fe: 



2/xr n r n _ 1 sin n sin 0„_ 1 ft 
Now by using the Taylor development up to order 4 



V (r n , 9 n ) e 



we have 

K(f b ,f a ;T) = 



cos A0 n = l~\ (A0„) 2 + ^ (A0 n ) 4 + •• 



1 



2nr b r a y/sm9 b sin 9 a 

N-l .N-l N 



/n*./ri*n(^)V^s) 

n=l n=l n=l 



A? 



AT-1 , AT +00 

n^nz «p* 

n=l n=l m=— 00 



mhAtpn + 



/i (Ar T , 



2 £ 



// (Afl w )' 



( 



m 2 — |) ft 2 £ 



24 £ 



2 / ur n r n _i sin n sin n _i 



- V (r„, n ) £ 
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According to McLaughlin- Shulman procedure [27], the term with (A^) 4 leads to a quantum 
correction that can be calculated with the help of the property 



N-l r N 

l[d6 n = / Y[d(A6. 

n=l J n=2 



(22) 



and the use of the integral 

du u 4 exp [zau 2 ] = I 



(- 

\4a 2 



du exp [iait 2 ] 



This correction is 



ihe 



Then the propagator K (f b ,f a ; T) takes the following path integral representation 
K(f b ,f a ;T) 



(23) 



(24) 



2nr b r a a/ sin 9 b sin 9 a 

N-l r N-l N 

n /' /n^ir'' 



n=l n=l n=l 

N-l , N +oo 

n=l n=l m=—oo 



2m he 



2inhe 



N 



u(Ar^f u (A6 n y 
2 e 2 e 



h 2 



(m 2 — |) h 2 e 



o ij ■ n --V{r n ,9 n )£ 

2ur n r n _ 1 sm9 n sin n _i 



(25) 



At this level, we see that only the integration over (p is straightforward. One cannot do 
integration over 9 and r because of the position dependent kinetic term ^r n r n _i (A9 n ) 2 and 
the term h 2 (2ur n r n _ 1 sin 9 n sin 9 n _i)~ 1 . Let us in that case do integration over ip. We obtain 



K(f b ,f a ,T) = 



+°° -im(tpt,—<p a ) 



2tt 



-K m (r b , 9 b ; r a , 9 a ; T) 



(26) 



m=— oo 



where the novel propagator K m (r b , 9 b ;r a ,9 a ; T) is given by 



K m (r b , 9 b ; r a , 9 a ; T) 



1 



r b r ay /sm6 b sin 9 a 



N-l 



n=l 



N-l 



N 



\ i <!<■„ / n^n 



n=l n=l 



2inhe 



—) 

2mhe / 

2 t2 



A' 



(Ar n ) 2 /x (A^ n ) z , & 2 

1 — r n r n _! 1 e 

2 e 2 e 8/zr n r n _i 



(m 2 — |) ft 2 £ 



2\xr n r n _ x sin 6> n sin 



- V (r n , 6> n ) e 



(27) 
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Here we note that the terms z — and 7. ^ — =—5 — represent an effective potential 

8/w„r„_i 2/ir n r n _i smB n sin0 n _i ^ r 

that describes the quantum corrections resulting from the passage to spherical coordinates. 
In the next step we separate the angle dependence from the radial one by using a simple 
time transformation. 



III. SEPARATION OF VARIABLES 

Having formulated the problem of angle-dependant anharmonic oscillator in the frame- 
work of path integrals in spherical coordinates, let us proceed to find exact solutions by doing 
separation of variables. First we define the fixed energy amplitude K m (r b , 9 b \ r a , 9 a ; E) to 
be the Fourier transform of K m (r b , 9 b ; r a , 9 a ; T) 



with 



K m (r b , 9 b ; r a , 9 a ; T) 



-00 

—e^ T K m (r b , 9 b ; r a , 9 a ; E) 



(28) 



K m (r b , 9 b ; r a , 9 a ; E) 



1 



r b r a y/ 'sin 9 b sin 9 a 



+00 

dTe l f 
00 ■' ,<-i 



exp 



Hdr n I n^n 

n=l n=l n=l 

Ei A* ( Ar n) 2 , H (A6> n ) 2 
I 7T 1" 7T r " r n-l h 



Imhe 



2 e 



h 2 



(m 2 — |) Tie 



2/xr n r n _i sin6> n sin 



- V (r n , 6> n ) e 



(29) 



In order to be able to separate angular motion from the radial one we use the method of 
[28] . To begin we change the evolution time from t to s, with 



ds 



dt 



(30) 



This transformation is equivalent to 



e (^n^nTn—l 



(31) 



with a time interval 



S 



dt 



(32) 



To incorporate this changes in the path integral representation of K m (r b ,9 b ;r a ,9 a ; E) , we 
start form the identity 

dSd ( S I ^ ] .1 (33) 



and by the use of 



(34) 
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we get 



dSr 2 (T)5[T 



r 2 ds = 1. 



(35) 



By inserting the later identity in (J29J) and being aware of r (T) = r&, we can write 
K m (r b , 9 b ; r a , 9 a ; E) as an integral of two independent kernels 



with 



+oo 



K m (r b , 9 b ; r a , 9 a ; E) = dS P E (r b , r a ; S) Q m (9 b , 9 a ; S) 



(36) 



where 



and 



Pe (r b ,r a ; S) = 

N-l N 

n^n 



n 

r„ 



n=l n=l 



2mho n r n r n _ 1 



exp 



E 



\i (Ar>, 



2 o n r n T n — i 



-mu 2 r 2 n + 



(a - p) re 

2/< 



r n r n _icr n + (E + V ) 



8/i 



Or, 



(37) 



1 



V sin 6*;, sin 9 a 



exp 



E 



JV-l N 

n^n 

n=l n=l 
2 



2inha n 



n(M n y h 2 f\ 2 -{ + k 2 - 



2 <y n 2/i \ sin 2 # n cos 2 9 n 



Or, 



(38) 



k = \H + 



(39) 
(40) 



In the continuous limit the angular propagator Q m (9 b , 9 a ; S) has the form 



Q m (e b ,0a,s) = -=l== [do exp'- [ S (^0 2 -f (^f + ^i))^ (4i) 

\/sm9 b sm9 a J nj \2 2/^ \ sin 9 cos 2 9 J J 

which has some resemblance to the Poschl- Teller problem. Q m (9 b ,9 a ; S) is then integrable 
and the result is I29H3T1 



Qm (9 b , 9 a , S)=J2 e " H " 9)5 e no,m (*a) 6„ fl , m (9 b ) (42) 
n e =0 



where e (n e ) is given by 



h 2 

e(n e ) = —(2n e + k + X + l) 2 (43) 
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and angular wave functions are 

®n e , m (e)=N„ (sin#) A (cosfl)^ P^(cos20) 
with the normalization constant 



N, 



ng,m 



l 2(2n e + k + A + l) 



n e \T (n e + k + A + 1) 
r (n g + k + 1) T (n e + A + 1) 



(44) 



(45) 



For P E (r b ,r a ; S) we can't do integration directly because the kinetic term „— — ( ' A J' n ' > 
contains an inconvenient variable dependent mass. It is the necessary to restore the original 
time t. To this aim let us proceed as follows. First we incorporate the result of Q m (8 b , a ; S) 



in (36) to get 



K m fa, 6 b ; r a , 6 a ] E) = ^ Q n g ,m (°a) @n e ,m fa) ^n g ,m fa, T a ; E) 

n e =0 



where K, ne , m fa, r a ; E) is given by 



/•+oo 

iC ne , m fa, r a , E)= dS e~^ s P E fa, r a , S) 
Jo 



and has the following path integral representation 

N-l N 



K"n,g,m (j'bi 1 a, Ej 



n 



r. 



a JO 



dS / 11^,, If 



exp S 



ji (Ar. 



n=l n=l 
2 



2 @rJ'rJ'n—\ 



2mha n r n r n _i 
1 2 2 , {a-P)h 2 

2 mu r « + 



VnTn—l&n 



h 2 



+ {E + Vq) r n r n ^a n - ( %- (2n e + k + A + l) 2 - — i a„ 



(46) 



(47) 



(48) 



Then we make the change r n r n ^\o n — > £ to obtain for 1C ne , m {r b , r a ; E) the following form 



1Cng,m fa,r a ; E) 



+oo 



. d T / n^n 

1 b 1 a Jn J -, , 

n=l n=l 

x2 



AT-1 AT 



exp 



h 



/x (Ar„ 



2 £ 



2m he 



, o (a-/3)h 2 



2 mUrn+ 2^rl 

2 1 



(£ + V )e- 



fi 2 (2n e + fc + A + 1) - 
2^ ^~ 



(49) 
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For K, ng , m (r b , r a ; T) we have 



!Cn e ,m (n, r a ;T) = — / Yl dr n JJ ( 

a n=l n=l ^ 



2mhe 



exp 



E 



A* (Ar n )' 



\ 2 e 



-V + -muj 2 r 2 + 



(a - p) h 2 



2 1 



h 2 (2n e + fc + A + 1) — 
2/1 



(50) 



In the next section we show how to find solutions by integrating /C ne/m (r&, r a ; T) . 

IV. INTEGRATION OF THE RADIAL PROPAGATOR 

Having shown how to do separation of variables and integration over angular ones let us 
do integration over radial variable to find the final solution to our problem. Starting from 
the K, ng ,m (?"&, r a ] T) that can be written in the form 



1 



fcn e (n, r a ] T) = — ^— / ] [ dr„ J [ 



N-l 



N 



n=l n=l 
2 



2mhe 



i \ / n (Ar, 



n / 



e 



h 2 

-V5, + + - 




(51) 



with 



/ = V /(A; + A + 2n e + l) 2 + (a - (3) - i 



(52) 



This path integral is the ordinary path integral of the radial oscillator, the solution is known. 
We can find in [B] that K ne (r b , r a ; T) takes the form 



fcn g (r b ,r a ;T) 



iV T ^ u 



sjr h r a ihsin(uT) 



\ hi sin (uT) 



exp 



(53) 



Now to obtain the radial wave functions and the energetic levels we use the Hille-Hardy 
formula 



l-s : 



exp 



2\[XYi 
1 - s 2 



°° 2n+£+§ TJ | (n -i(X+Y) , , xJ 



ra=0 



T (n + a + 1} 



"2 



(54) 
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and we take s = e iUjT ; X = ^r\\ Y = ^fV^, to get the spectral decomposition of 

K-"ng,m (j'b; '"aj -^J 



n=0 



-i[u{2n+l+^ )-Vo\T_ 



TV. 



r[n + £ + 



or in closed form 



JC ng ,m {r b , r a , T) = ^ e n E n,n e , m T x ^ ^ 



(56) 



n=0 



where the energy spectrum is given by 



E 



n,ng,m 



2n + t + - ) hw - V 



(57) 



and the radial wave function reads 
R 



<<n,ng,m 



(r) 



UU 2 

exp,--r 



\ V h ' r(n + £+fj 
Finally the relative propagator can be expressed in the following spectral decomposition 

oo oo +oo 

k (f b , f a , t) = J2J2 E e ~ iEn ' ne ' mT r n , ne ,m (fa) * 

where 



n=0 ng=0 m=—oo 



J n,ng,m 



in 



(59) 



ipnim (r) = N ntns>m exp ( -T^r 2 




x (sin 0) X (cos d) k+1 * (cos 29) e 

with the normalization constant 



limp 



N, 



n,ng,m 



\ 



2 ///OA 



7T 



V h ) 



(2n e + k + A + V 



n\n e \T (rig + k + A + 1) 



r (ng + k + 1) T (n* + A + 1) T In + I + 



(60) 



(61) 



Let us remark that when we consider the particular cases studied in [T9TI2T] , we see that our 
results coincide with these results. 



V. CONCLUSION 

In this paper we have given a straightforward method to solve the problem of noncentral 
anharmonic oscillator in three dimensions. In the first stage we have expressed the relative 
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propagator by means of path integrals in spherical coordinates. Then by making an adequate 
change of time we were able to separate the angular motion from the radial one. The angular 
part path integration is reduced to the well-known Poschl-Teller problem and the resulting 
radial path integral is written in the form of three dimensional isotropic harmonic oscillator. 
Then we have exactly calculated the relative propagator and we have extracted the bound 
states energies and the corresponding wave functions. 

We remark also that the presented method is useful for all potentials of the form 

V(r,8)=v(r)+ ^ (62) 

and there is no need to use other systems of coordinates. 

Through the formulation given above and the obtained energies and wave functions we 
conclude that the path integral formulation is a powerful method to study quantum dynamics 
of particles in nonrelativistic theory. 
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